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Because of the inherent nonlinearity, most fluid-dynamical problems must be solved by ei- 
ther analytical approximations or by numerical computations. In this chapter we shall first 
explain the lubrication approximation, with applications to the flow of thin layers. Atten- 
tion will then be turned to the slow flow past a sphere and a cylinder. Application to the 
environmental pblem of aerols will be briefly discussed. Lastly the slow withdrawl of fluid 
from a straitified reservoir into a sink will be analyzed. 



2.1 A thin fluid layer flowing down an incline 

As our first example, we consider the flow of a thin layer of viscous fluid on an inclined plane. 
Referring to Figure 2.1.1, let the coincide with the plane bed inclined at the angle 

9 with respect to the horizon, and the y— axis normal to the plane bed, 




Figure 2.1.1: A thin fluid layer flowing down an incline plane 



2.1.1 Exact equations 

Continuity : 



Momentum 



and 



du du du 
dt dx dy 



dv dv dv 
dt dx dy 



du dv g 
dx dy 
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= g sin 9 - -— + v — + — 
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(2.1.1) 
(2.1.2) 
(2.1.3) 
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2.1.2 Exact boundary conditions 

On the rigid bottom no slipping is allowed, 

u = v = 0, y = 0. (2.1.4) 

The top of the fluid layer of local depth h(x, t) is a material surface whose position is un- 
known until the problem is solved. It is necessary to specify a kinematic boundary condition 
concerning the motion, and two dynamic conditions concerning the stresses. 
Kinematic condtion: 

Since the top is a material surface, it must be always composed of the same fluid particles. 
The fluid velocity must be tangential to the moving surface at all times. Let us describe the 
top surface by 

F(x, t) = F(x, y,t) = y- h(x, t) = (2.1.5) 

Consider a surface point which is at x at time t, then it will move to x + q s dt at time t + dt 
where dt is a very short time interval and q s is the velocity of the surface point. The new 
top surface is then described by 

F(x + q s dt, t + dt) = 
Taylor-expanding for small dt, we get 

F(x, t) + dt + Q, ■ VF(f , t)) + O(dt) 2 = 

The first term vanishes because of (2.1.5), hence 

^ + a-vF(f, t ) = o 

in the limit of dt — 0. Now the unit normal to the surface is 

VF(x , t) 



n 



|VF(£,t)| 



and the normal component of the surface velocity is q s ■ n. Since the normal velocity of a 
fluid particle q on the surface must be equal to the normal velocity of the surface at the 
same point, 

g s ■ VF(x,t) _ q- VF(x,t) 
\VF(x,t)\ ~ \VF(x, t)\ 

it follows that 

^^ + g-VF(f,t) = (2.1.6) 

This is the kinematic bounday conditon on the top fluid surface. Substituting F = y — 
h(x, i) = we get the kinematic boundary condition. 

dh dh , . . _ „ „. 

dt +u dx~ = v > y = h ^ (2 - L7) 
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Dynamic boundary conditions 
Let the free surface be stress-free: 



'-pfiij + nj) rij = 0, i = 1,2. (2.1.8) 



or 



51 = —pn-i + Tnni + T U n 2 = (2.1.9) 

5 2 = -pn 2 + T 2 in 1 + r 2 2n 2 = 0. (2.1.10) 
on y = h. The unit normal vector has the components: 

dh i 

"i = - , te . 2 , "2= , = (2.1.11) 



So far these exact boundary conditions are highly nonlinear, in part because the position 
for applying these conditions are not known. As an useful consequence of (2.1.7) to be used 
later, we integrate the continuity equation in depth, to get, after invoking Leibnitz's rule 

J \'q^ + 'q^)^ z = 'q^j u dx-u(x,h,t)— + v(x,h(x,t),t) = 

With the help of (2.1.7), we get the depth-integrated law of mass conservation, 

dh d ( fh \ dh dQ dh d(uh) n , ft1in , 

where Q denotes the rate of volume discharge and u the depth averaged horizontal velocity 

B =£ = 5jf»* (2 ' Ll3) 



